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I.  INTRODUCTION 


In  1970,  Ward  and  Mansfield1  made  some  hypersonic  damping- in-pitch 
measurements  of  9-degree  cones  by  the  free  oscillation  technique.  Since 
an  air  bearing  support  was  used,  the  models  were  free  to  spin  as  well  as 
pitch  and  yaw.  An  apparently  minor  tip  asymmetry  was  induced  by  a 
planar  cut  40°  off  the  X-Z  plane.  The  projected  frontal  area  of  the  cut 
was  approximately  0.5%  of  the  base  area. 

The  free  oscillations  were  then  fitted  by  the  Nicolaides  tricyclic 
theory.2  This  theory  assumes  that  the  sole  effect  of  small  asymmetry  is 
to  add  a missile-fixed  constant-amplitude  trim  moment  term  to  the  moment 
expression  for  the  unmodified  symmetric  cone.  A missile  that  satisfies 
this  assumption  is  said  to  have  a slight  configurational  asymmetry.  It 
was  found,  however,  that  the  measured  motion  cannot  be  reasonably  rep- 
resented by  the  tricyclic  theory.  In  1972,  Walchner,  Sawyer,  and 
Yelmgren3  showed  that  the  apparently  very  small  asymmetry  used  in  the 
AEDC  tests  was  sufficient  to  make  the  pitch  and  yaw  frequencies  unequal 
and  thus  a more  complex  mathematical  moment  assumption  is  required  to 
analyze  the  free  oscillation  data.  References  4 and  5 discuss  two 
other  cases  of  apparently  symmetric  missiles  whose  motions  seem  to  be 
poorly  described  by  either  the  epicyclic  motion  performed  by  symmetric 
missiles  or  the  tricyclic  motion  performed  by  slightly  asymmetric 
missiles. 


2.  L.K.  Ward  and  A.C.  Mansfield , "Dynamic  Characteristics  of  a 9- 
Degree  Cone  with  and  without  Asymmetries  at  Mach  Number  10 , " 
Arnold  Engineering  Development  Center  TR-70-1,  March  1970. 

2.  J.D.  Nicolaides,  "On  the  Free  Flight  Motion  of  Missiles  Having 
Slight  Configurational  Asymmetries ,"  Ballistic  Research  Labora- 
tories Report  No.  858,  June  1955,  AD  26405.  (See  also  IAS  Pre- 
print 395,  1952.  ) 

3.  0.  Walchner,  F.M.  Sawyer  and  K.E.  Yelmgren,  "Aerodynamics  of  a 
Slender  Cone  with  Asymmetric  Nose  Bluntness  at  Mach  14,"  AIAA 
Journal  10,  August  1972,  pp.  1121-1122. 

4.  P.F.  Intrieri , D.B.  Kirk,  G.T.  Chapman  and  J.E.  Terry , "Ballistic 
Range  Tests  of  Ablating  and  Nonablating  Slender  Cones,"  AIAA 
Journal  £,  March  1970 , pp.  558-564. 

5.  P.  Jaffe,  "Nonplanar  Tests  Using  Wind  Tunnel  Free  Flight 
Technique ,”  Journal  of  Spacecraft  and  Rockets  10,  July  1973, 
pp.  435-442. 


The  effect  of  constant  axial  spin  on  the  stability  of  aircraft  was 
first  considered  by  Phillips.6  Recently,  Hodapp7-8  has  extended  this 
work  to  unsymmetric  re-entry  vehicles  in  a most  elegant  fashion.  Both 
Phillips  and  Hodapp  make  use  of  missile-fixed  coordinates.  These  co- 
ordinates are  convenient  if  the  missile  is  not  spinning  or  if  on-board 
instrumentation  is  used.  For  the  three  experiments  cited  earlier,  ex- 
ternal photographic  instrumentation  was  used  and  the  nonrolling  aero- 
ballistic  axes  are  much  more  convenient.  A second  difficulty  of  the 
Phillips  and  Hodapp  papers  is  that  the  general  pitching  and  yawing 
motion  of  unsymmetric  missiles  is  described  by  five  modal  vectors  and 
the  results  are  given  in  a form  in  which  it  is  difficult  to  see  how 
the  five  vectors  reduce  to  the  three  modal  vectors  appropriate  to 
symmetric  motion  with  trim. 

In  this  paper,  we  make  use  of  the  complex  variables  that  consider- 
ably simplify  the  analysis  of  symmetric  missile  motion.9  These  vari- 
ables allow  us  to  move  easily  from  missile-fixed  axes  to  nonspinning 
axes.  They  also  allow  us  to  obtain  very  simple  expressions  for  an 
"almost"  symmetric  missile.  All  results  are  obtained  for  constant  spin. 

In  Sections  II  through  VI,  the  analysis  is  limited  to  a linear 
yaw  moment  and  a linear  pitch  moment  with  C not  equal  to  - C 

ol  B 

Since  the  pitch  and  the  yaw  frequencies  are  different  for  zero  spin, 
there  is  a resonance  region  for  spin  when  the  spin  is  between  these 
frequencies.  After  deriving  exact  relations  for  the  two  frequencies 
of  the  combined  pitching  and  yawing  motion  for  a spinning  missile,  we 
obtain  approximations  for  an  "almost  symmetric"  missile.  An  almost 
symmetric  missile  is  one  whose  zero-spin  pitch  frequency  is  very  near 
its  zero-spin  yaw  frequency.  For  these  approximations,  three  regions 
of  spin  are  considered  separately:  spin  near  zero,  spin  near  reso- 
nance, and  spin  far  from  zero  spin  and  resonance  spin.  For  the  third 
region,  it  is  shown  that  in  the  aeroballistic  system,  two  of  the 
frequencies  are  approximated  by  the  fast  and  slow  rates  <^iS  , 4>2S 

appropriate  to  a symmetric  missile  with  a static  moment  derivative 


6.  W.H.  Phillips , "Effect  of  Steady  Rolling  on  Longitudinal  and 
Directional  Stability,"  NACA  Technical  Note  No.  1627,  June  1948. 

7.  A.E.  Hodapp,  Jr.,  "Effects  of  Unsymmetrical  Stability  Derivative 
Characteristics  on  Re-Entry  Vehicle  Trim  Angle  Behavior, " 

Journal  of  Spacecraft  and  Rockets  11,  May  1974,  pp.  300-307. 

8.  A.E.  Hodapp,  Jr.,  "Effects  of  Unsymmetrical  Stability  Derivative 
Characteristics  on  Re-Entry  Vehicle  Transient  Angular  Motion," 
Journal  of  Spacecraft  and  Rockets  13,  February  1976,  pp.  82-90. 

9.  C.H.  Murphy,  "Free  Flight  Motion  of  Symmetric  Missiles ,"  Ballistic 
Research  Laboratories  Report  No.  1216,  July  1963,  AD  442757. 
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(1.1) 


- Cn  )/2 


Two  other  modes  exist  with  frequencies  2p  - 4^  , 2p  - <P2$  . The 
amplitude  of  these  modes  is  shown  to  approach  zero  as  C -*•  -C 

a e 


Near  resonance  the  fast  mode  and  its  associated  frequency  2p  - $ 

both  approach  the  spin  rate,  p.  For  spin  rates  in  the  resonance  zone 
between  the  pitch  and  yaw  frequencies,  both  ^ and  2p  - 4^  become  equal 

to  the  spin  rate  and  one  mode  is  exponentially  damped  while  the  other 
is  exponentially  undamped. 

In  Section  VII,  the  effect  of  linear  aerodynamic  damping  is 
studied.  Its  major  effect  is  to  reduce  the  region  of  spin  rates  for 
which  exponential  undamping  exists.  Outside  the  resonance  region,  the 
two  damping  rates  are  very  similar  to  those  for  a completely  symmetric 
missile. 


II.  EQUATIONS  OF  MOTION 

The  velocity  equation  for  nearly  horizontal  flight  can  be  easily 
written  in  terms  of  the  aerodynamic  force  along  the  trajectory--the 
drag  force: 

mV  = - (PV2/2)SCD  (2.1) 

If  the  independent  variable  is  changed  from  time  to  dimensionless  arc- 
length,  s,  this  reduces  to 


The  usual  coordinates  to  describe  the  angular  motion  of  a rigid 
body  are  body-fixed  coordinates.  For  the  case  of  a symmetric  missile, 
the  aeroballistic  axes,  which  pitch  and  yaw  with  the  missile  but  have 
a zero  spin  rate,  are  more  convenient.  In  this  study  of  an  "almost" 
symmetric  missile,  we  will  use  missile-fixed  axes  to  obtain  our 
results  and  then  transform  the  results  to  the  aeroballistic  axes. 

The  X-axis  is  fixed  in  the  body  along  the  principal  axis  of 
inertia  nearest  to  the  desired  flight  direction.  Since  we  retain  the 
assumption  of  mass  rotational  symmetry,  the  transverse  moments  of 
inertia  are  equal  and  for  any  roll  orientation  of  the  Y and  Z axes, 
these  axes  will  also  be  along  principal  axes  of  inertia.  The  actual 
orientation  of  these  axes  will  be  selected  to  simplify  the  aerodynamic 
moment.  (The  effect  of  small  mass  asymmetry  is  given  in  Appendix  A.) 

11 


The  equations  of  motion  can  be  written  in  terms  of  the  aerodynamic 
force  and  moment  and  the  derivatives  of  the  linear  and  angular  momentum. 
(For  simplicity  the  effect  of  gravity  has  been  neglected.  Its  effect 
on  the  angular  motion  is  usually  quite  small.) 


m V = (p V2/2) S (Cx,  Cy>  C2) 


fi  = (pV2/2)S*.  (Cz,  Cn,  Cn) 


where 


V = (u,  v,  w) 


» = ux  P.  It  q>  It  r) 


The  vector  derivatives  can  be  computed  by  the  well-known  relation  for 
any  vector  (U^,  Uy,  Uz) : 

(UxJ  V “z3’  = V + “ * (Ux'  Uy>  Uz)  (2‘5) 


where 


“ = (p,  q»  r) 


Equations  (2. 3-2. 4)  can  be  expanded  by  use  of  Equation  (2.5)  and 
the  independent  variable  changed  from  t to  s.  For  small  angles  (u  ~ V) , 
the  transverse  components  of  these  vector  equations  become 


v'  . rl  , f w\  _ , 

~ + v~  ' * { v)  ' S 


^ * (i  - ^ (f)  - < 
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The  complex  notation  that  is  so  helpful  for  a symmetric  missile  can  now 
be  introduced.  The  complex  angle  of  attack  and  transverse  angular 
velocity  have  the  usual  definitions: 

5 = ~Y--W-  = 6 + i a (2.10) 

p = (q  + i r)  l/V  (2.11) 


By  differentiating  Equations  (2.10)  and  (2.11)  we  obtain  (with  the 
aid  of  Equations  (2.2)  and  (2. 6-2. 9)): 


C'  + i - i p 


Ccv 


i C*) 


+ CD? 


(2.12) 


p'  + i (<j>' 


P)p 


«£  + 1 <£>  + CD  * 


(2.13) 


where 


p - ax/it)*' 

k2  = It/m^2 


III.  LINEAR  STATIC  MOMENT 

The  most  general  moment  components  that  are  linear  in  a and  B 

are 


C = C + C a + C B 
m m„  m m„ 

0 a B 


(3.1) 


C 

n 


C + 
n 

0 


C a + 
n 


(3.2) 
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In  complex  form,  this  becomes 


1 


l 

i 


C + i 
m 


C 

n 


C + 
mo 


e 


+ 


(3.3) 


where 


The  effect  of  rotating  the  Y and  Z axes  through  the  angle  60  can  be 
easily  calculated  by  multiplying  Equation  (3.3)  by  exp  (i 0 Q)  . 


! 

c 

i 

i 


(3.4) 
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According  to  Equation  (3.4),  the  coefficient  of  the  complex  angle  is 
unaffected  by  a rotation  of  coordinates  but  the  coefficient  of  the 
conjugate  of  the  complex  angle  is  multiplied  by  exp  (2i60).  This 
property  allows  us  to  choose  a convenient  orientation  of  the  Y and  Z 
axes.  This  orientation  will  be  selected  so  that  the  coefficient  of  the 
conjugate  complex  angle  is  zero  or  a negative  imaginary  number 
(Cms  = 0,  < 0).  In  terms  of  the  coefficients  defined  in  Equations 

a a 

(3. 1-3. 2), 


C = C 
m0  n 

6 a 


-C  > c 
m = nfi 
a p 


For  a symmetric  missile,  C,.  is  zero  and  C„„  is  an  odd  function 

M Mb 

a a 

of  spin.  (The  side  moment  for  this  case  is  usually  called  a Magnus 
moment.)  The  addition  of  a complex  constant  (C  + i C ) to  the 

p m0  n0 

moment  of  a symmetric  missile  introduces  a trim  angle  and  the  possibility 
of  resonance.2'  10-11 

The  primary  effect  of  the  CMg  side  moment  coefficient  on  the 

a 

angular  motion  of  a symmetric  missile  is  on  the  damping  rates  and  not 
on  the  frequencies.  The  damping  rates  are  assumed  to  be  small  compared 
with  the  frequencies  and  will  be  calculated  later  as  a perturbation  of 
the  zero-damped  motion.  Thus,  the  side  moment  will  be  grouped  with  the 
aerodynamic  forces  and  the  damping  moments  and  neglected  until  the 
perturbations  are  introduced.  Hence  the  static  moment  to  be  used  in 
the  first  stage  of  our  analysis  is  defined  by 


C + i C = C 
m n m„ 


1 Cn  ' 1 fCM  5 ' 

0 La  a 

+ C a + i|C  + C 
ma  |_  n0  n£ 


10.  R.  L.  nelson,  "The  Motions  of  Bolling  Symmetrical  Missiles 
Referred  to  Body-Axis  System ,"  NACA  Technical  Note  No.  1627, 
November  1956. 


11.  D.A.  Price,  Jr.,  and  L.  E.  Ericsson,  "A  New  Treatment  of  Roll- 
Pitch  Coupling  for  Ballistic  Re-Entry  Vehicles AIAA  Journal  8 , 
September  1970,  pp.  1608-1615. 


IV.  FREQUENCIES  AND  MODAL  AMPLITUDES 

The  static  moment  of  Equation  (3.5)  can  now  be  inserted  in 
Equation  (2.13),  the  aerodynamic  force  neglected  (C^  = Cy  = = 0) 

and  Equation  (2.12)  used  to  eliminate  y . We  obtain 


where 


5"  + i f 4>'  s'  ■ C + M ^ = i A 


I 

f = 2 - -p 


M = M + (f  - 1)(*')2 
r 


(4.1) 


M 


M 


eM  < 0 
a 


Since  Equation  (4.1)  is  a fourth-order  differential  system  in  the  real 
variables  a and  6,  periodic  solutions  correspond  to  two  pairs  of  conju- 
gate roots  of  a fourth-degree  polynomial. 

For  a symmetric  missile  with  trim,  M is  zero  and  the  solution  to 
Equation  (4.1)  can  be  written  in  terms  of  two  rotating  complex  vectors 
plus  a constant  trim  vector. 


C = 


iu>  s 
e 1 


+ k 

2 


ioj  s 
e 2 + 


k 

3 


(4.2) 


where 


U)  . 


J 


[(f  *')2 


[P2 


k - - i A/m 
3 r 


In  the  nonrolling  aeroballistic  coordinates,  this  becomes 


£ = £ e 


i<j/  s 


= kj  e 


i<t>'s  . 

1 + k e 

2 


i<t»'s  , 

2 + k3  e 


• J 
i<t>  s 


(4.3) 


where 


4/.  = u . 
J J 


+ 


M>  ± [P2  - 4M] 


Note  that  the  trim  vector  becomes  infinite  for  «r  - 0 . This  is  the 

usual  spin-yaw  resonance  condition  and  the  spin  that  makes  M vanish  is 
resonance  spin. 

Usually  <j/  is  defined  to  be  the  larger  frequency  and,  if  we  limit 
the  spin  to  positive  values,  Wj  would  then  correspond  to  the  positive 
square  root  and  to  the  negative  square  root.  This  means  that  for 
negative  M,  u is  always  negative  and  w is  positive  for  M negative 

2 lx 

and  negative  for  positive.  Thus,  u>1  changes  sign  at  resonant  spin 
and  is  always  less  in  magnitude  than  u . 


I 


Since  Equation  (4.1)  contains  £,  the  solution  should  consist  of 
complex  vectors  and  their  conjugates. 


£ = kx 


iw,  s 

e 1 + k2  e 


+ k 


3 


-iu  s -iw  s 

+ ki,  e 1 + k5  e 2 

By  direct  substitution  into  Equation  (4.1),  we  see  that  the  trim 
vector  takes  the  very  simple  form 


(4.4) 
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(4.5) 


- i[Mr  A - M A] 

M2  - M2 
r 


Resonance  now  occurs  for  the  two  values  of  spin  that  make  the  denomi- 
nator in  Equation  (4.5)  vanish.  In  terms  of  pitch  and  yaw  zero-spin- 
rate  frequencies  and  u»  , these  spin  rates  are 


♦i  - ^ 


i = a,  6 (4.6) 


where 


o.  = [ - k”2  C*  ]“* 
a t m 

a 


r.  “2  r* 

= [kt  Cn 

p 


Note  that  the  pitch  and  yaw  axes  were  selected  so  that 


t'  > 


(4.7) 


Direct  substitution  of  Equation  (4.4)  into  Equation  (4.1)  yields 
two  pairs  of  equations,  for  k and  k^  and  for  k2  and  k5  . The  first 

pair  can  be  written  in  the  form 


where 


k [C  + 2 f <f>'  iu1  ] - k4  M = 0 

-kM  + k C =0 
i 4 1 

C = u)2  - f <j/  a)  + M 
1 1 1 r 


The  solution  to  Equations  (4. 8-4.9)  is 


(4.8) 

(4.9) 

(4.10) 


(4.11) 
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(4.12) 


The  frequency  equations  follow  from  Equations  (4.10,  4.12)  and  the  ' 

corresponding  equations  for  the  second  frequency.  j 

i 

I 

u1!  - 2 a si?  + b = 0 (4.15)  • 

J J I 

j = 1,  2 ] 
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where 


a = (1/2)  (f  <j/)2  - M 


b = M2  - M2 
r 


= (f  - 1)2[(*')2 


(<)2][(^)2  - 


The  solution  to  Equation  (4.15)  can  be  written  in  the  form: 


2 

cor  = a 


± s/a 


J * 


- vT 


(4.16) 


Since  a is  always  posit:  e,  Equation  (4.16)  gives  positive  values  of 

w2  and  a)2  for  b and  (a2  -b)  positive.  For  negative  b,  only  u2  is 
1 2 2 
positive  and  a real  does  not  exist.  Under  these  conditions,  the 

precise  relations  for  real  &>.  are: 

" ' vT^ 


u =6 
1 


00  - 
2 


- fa*-)" 


(4.17) 


(4.18) 


If  Equation  (4.11)  and  its  k5  counterpart  are  substituted  in 
Equation  (4.4)  and  if  the  kjs  are  expressed  in  the  polar  form 
Kj  exp  (i  ^e  solution  motion  takes  a very  simple  form  in  aero- 

ballistic  coordinates: 


r i<J>.  - -i  i (2<f>  - <f>j)  . 

C = [e  1 + M Cj  e ] 

1 i(24>  - 4>2) 


K2  [e  2 + M C e 


(4.19) 


+ k eU 

3 
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i 


r 


for  0 < 4/  < <(>g  or  4/  > 4^  , where 

4>  =4,  + <t>  ! S 

9j  93  0 93 

4(  “ u) . + 

3 3 

4>  = (J/S 


3 = 1.  2 


V.  RESONANCE  REGION 


When  the  spin  rate  4/  is  between  4>g  and  4>^  , b is  negative  and 

Equation  (4.17)  is  invalid.  (It  should  be  noted  that  Equation  (4.18) 
does  give  a real  number  for  ui  although  the  calculation  involves 

complex  numbers.)  For  this  resonance  region,  the  terms  with  u]  in 
Equation  (4.4)  must  be  modified: 


where 


r . X s . -Xs  . ^U)2S 

5 = klRe  + k4Re  ♦ k2  e 

-iu>  s 

+ k5  e ^ + ^ 

X > 0 


(5.1) 


k.R  = K.r  exp  (i  4jR),  j - 1.  4 

Substitution  of  Equation  (5.1)  in  Equation  (4.1)  yields  two  pairs  of 
real  equations  for  the  coefficients  of  exp  (Xs)  and  exp  (-Xs). 


(X2 

- M 

r 

+ ft) 

cos 

4 - f 4/ 

iR 

X sin  4 B 

I K 

0 

(5.2) 

f 

4/  X 

cos 

$lR  + 

ho 

1 

n2 

1 

M)  sin  $lR  = 

0 

(5.3) 

(X2 

- M 

r 

+ M) 

cos 

4>  + f 4/ 

4R 

X sin  4 s 

4 R 

0 

(5.4) 

f 

4'  X 

cos 

4 > + 

4R 

(X2  - M 

r 

- M)  sin  4 = 

4R 

0 

(5.5) 

21 

According  to  these  equations 


tan  4>  = - tan  4>, 

lR  4 R 


-f  X 


X2  - M - M 


X4  + 2a  X2  + b = 0 


The  real  positive  solution  to  Equation  (5.7)  is 


X = [ - a + (a2  - b)"5]"4 


For  an  almost  symmetric  missile  (M  <<  M) , b <<  a2  and 


a = (f  4>'  ) 2/ 2 


X = [M2  - Mz]h/(f  *')  = V-b/  (f  *') 


tan  4> 


(*')2  - C^)2  ** 
C^)2  - (*')2 


(5.10) 


(5.11) 


According  to  these  approximations,  X vanishes  for  spin  on  the  boundaries 
of  the  resonance  region  and  is  a maximum  near  the  center  (M  = 0) . The 

planes  of  the  exponentially  growing  and  decaying  angles  are  near  the 
pitch  plane  for  spin  near  <j>^  and  are  near  the  yaw  plane  for  spin  near 

. For  maximum  growth/decay  rates,  the  planes  of  the  growing  angle 

and  the  decaying  angle  are  perpendicular  to  each  other  and  both  planes 
bisect  the  angle  between  the  pitch  and  yaw  planes. 


r 


1 


VI.  C.  AND  w.  FOR  ALMOST  SYMMETRIC  MISSILES 
3 3 

We  will  now  obtain  special  relations  for  an  almost  symmetric 
missile  for  which  fil  is  much  smaller  than  M.  For  this  case,  the  radical 
in  Equations  (4.14)  for  can  be  expanded  for  three  regions  of  spin- 

near  zero,  near  resonance  and  otherwise.  These  expansions  yield  the 
results  in  Table  I. 

Far  from  resonance,  m£  >>  M2  and  a good  estimate  of  the  frequencies 
can  be  obtained  from  the  approximation 

VT  = 6 (Mr  - M2/2Mr)  (6.1) 


Equations  (4.17 — 4.18)  can  now  be  expanded  in  terms  of  M to  yield: 


fr'  = *'js  ♦ g.  M2  (6.2) 


= y [p  ± v/p2  - 4M  ] 


(6.3) 


where  the  g^'s  are  given  in  Table  II.  According  to  these  equations, 

very  good  approximations  to  the  frequencies  of  an  almost  symmetric 

missile  are  the  symmetric  frequencies  of  a missile  with  the  average 

moment  coefficient  Cw  of  the  almost  symmetric  missile.  This  is 

M 


a 

correct  to  within  a term  in  M2 


Note  that  the  coefficients 


83 


are 


unbounded  for  zero  spin  but  only  g^  has  difficulty  near  resonance. 
(g2  has  the  limit  * [(f  )"3]/2  as  Mr  -*■  0 .) 


Near  resonance,  Equation  (4.17)  can  be  expanded  for  small  b to 
give  the  relation  in  Table  II.  A little  additional  algebra  yields  the 
expansion  for  frequencies  near  zero  spin  given  in  Table  II. 

According  to  this  table,  § at  zero  spin  is  the  zero-spin  yaw 
frequency  and  <$/  increases  as  tiie  spin  increases.  At  frequencies  well 
away  from  zero  or  , <(/  is  well  approximated  by  <f/g  . In  the 

resonance  region,  it  is  the  spin  frequency  and  at  spins  well  above  the 
resonance  region,  it  is  once  again  well  approximated  by  the  symmetric 
missile  fast  frequency. 
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<(/  has  a much  simpler  behavior.  At  zero  spin,  it  is  the  negative 
2 

of  the  zero-spin  pitch  frequency  and  increases  as  spin  increases.  Well 

away  from  zero  spin,  it  is  well  approximated  by  4/ 

2 b 


In  Table  III,  the  various  approximations  for  C are  used  to 

construct  expressions  for  the  various  motions  for  an  almost  symmetric 
missile  with  no  trim  moment  (k3  = 0).  This  catalog  of  motions  starts 

at  zero  spin  with  the  expected  Lissajous  figure  formed  by  combining  the 
pitching  and  yawing  motion  with  different  frequencies  (Case  A).  For 
slow  spin,  we  have  the  more  complicated  form  of  Case  B. 


Case  C is  the  general  case  for  spin  well  away  from  zero  or 
resonance.  It  is  a four-mode  motion  with  two  modes  quite  similar  to 
those  associated  with  the  symmetric  missile  epicyclic  motion  with 
frequencies  4/  and  <f/  . The  two  additional  modes  have  much  smaller 

amplitudes  that  go  to  zero  as  the  missile  becomes  symmetric  (M  -+■  0) . 
Their  frequencies  have  the  rather  odd  form  2<t>/  - <p'.  . 


As  the  spin  approaches  the  resonance  region,  <f/  and  2<j/  - 4> 
become  equal  and  the  motion  associated  with  Kj  changes  character.  The 
amplitude  of  the  2 <j/  - 4/  mode  approaches  the  amplitude  of  the  <(/ 


mode  (Case  D) . At  the  singular  spin  value  of  ^ , a very  special 

secular  solution  is  found  (Case  E) . In  the  resonance  region  (Case  F) , 
the  two  modes  associated  with  <f/  both  rotate  at  the  spin  rate  but  one 

mode  grows  exponentially  and  the  other  decays  exponentially.  The 
exponent,  X,  depends  on  the  size  of  the  resonance  region  and  how  far 
the  spin  rate  is  from  the  boundaries.  Case  G is  the  secular  solution 
for  the  upper  boundary  of  the  resonance  region  and  Case  H is  the  motion 
for  slightly  higher  spins. 


In  order  to  show  this  spectrum  of  motions  in  a more  specific  way, 
the  frequencies,  modal  amplitudes  and  resonance  damping  rate  have  been 
computed  for  an  almost  symmetric  missile  with  C^  = 0.095  C^  and 

a a 

f = 1.9.  For  these  values. 


*'/*'  =1.10 

a B 

"«/*B  B 1-04 
wb/4/b  = 0,95 
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Table  III.  Aeroballistic  Motion 


lA.  . 0 


S ■ ( * H cos  9 ♦ 2 i K sin  9 

11  i 2 


■ *,[•“■  • 


B.  0 < $'  < e 


. i i(2*  - O 
il  - f ♦'  u M"  ] e J 


r i*2  , . ,,  i ^2*  - ♦,) 

* K2  [e  - [1  ♦ f ♦'  «2  M-l]-1  e 2 J 


A 1 <<  Af  A f 


1 • ” ” ’B*  ’ % 1 

M /->  c ./  ^-1  1 " Ol 

* M (2  f t wp  e 1 1 

r i*,  . -1  i(2,f  ■ *2^ 

♦ K2  e 2 - M (2  f 9'  u2)  e 2 J 

K 
^[e  * 


D.  - e < *'  < *' 


[1  - f ♦'  wj  M’1]-1  e 


i(2$  - ♦ 


‘1  * K’  [’* 


i* 


2 , -1  i(-2*  ' *J 

- M (2  f ♦'  u2)  e 2 


] 


5 = K e'^  I cos  9 * [i  ♦ 2 M (f  9' )' 1 s]  sin  9 
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E.  ♦'  = 


'i-l 


i r i4,  - \ i^2i>  ■ *2^  1 

J ♦ K le  2 - M (2  f ♦'  u2)  e j 


F.  < <fi  K (Resonance  Region) 


i*  T Xs  * i «lP  -As  - i*  D"| 

5"e  |>e  1R+K4Re  lRj 


[■ 


It2  , -1  i(2t  - ♦ ) 

e • M (2  f *'  u ) e 2 

2 


G.  i>/  - 91 
a 


5 » Kla  e1*  [1  + 2 i M (f  <,') 


i 1 h 

s]  cos  9 + i sin  * I + K le 

a aj  2 L 


i*.  - „ , -l  i(2<t  - O 

2 - M (2  f 9 e 2 


H.  9'  < 9'  < 9'  * z 


i*.  , ' 1 1 i(2*  - 4>,)1 

= Kj  e 1 - [1  + f 9'  Ul  M-1]  1 e 1 


r i*2  . 

* Kj  le  2 - M 


1 i (2*  - <t>  ) 
(2  f ♦'  u ) e 2 

2 


[T 

a'  >>  a' 

L_ 

a 

Same  as  Case  C 


where 


tan  «iR  = {[(♦') 2 - (»'8) 2 ]/[(♦;) 2 - (O')2]}’5 


and  where  K1I)(  KUn>  4 , 9.,  K,  , K are  parameters  determined  by  initial  conditions. 
ak  k a d la  ig 
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Figures  1 and  2 give  the  actual  frequencies  as  functions  of  4//4’/fi 

p } 

and  compare  them  with  the  corresponding  frequencies  for  a symmetric 

missile  with  moment  coefficient  C..  . 4/  is  essentially  4/_  except 

Ml  lb  ] 

oi 

near  zero  spin  and  resonance  and  is  essentially  <)/g  except  near  1 

zero  spin.  At  zero  spin,  4/  and  4/  reduce  to  to^  and  -to respectively.  | 

Figures  3 and  4 show  the  variation  of  the  "asymmetric"  modal 

amplitudes  K and  K in  comparison  with  the  corresponding  "symmetric"  j 

4 5 

modal  amplitudes  and  1<2  . K4  is  quite  small  except  near  zero 

spin  and  near  resonance  and  is  quite  small  except  near  zero  spin. 

It  is  important  to  note  that  is  always  a larger  fraction  of  Kj 

than  K is  of  K . This  is  due  to  the  presence  of  to.  in  C.  and  the 

5 2 1 J 

fact  that  u is  always  smaller  in  amplitude  than  to^  • Thus  we  would 

expect  to  observe  the  presence  of  the  24/  - 4/  mode  before  the 

2<)/  - 4>^  mode  has  a noticeable  effect. 

Finally,  the  resonance  damping  rate  is  plotted  in  Figure  5.  The 
peak  value  of  .047  corresponds  to  about  a 30%  increase  in  amplitude 
during  one  revolution. 
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VII.  EFFECT  OF  DAMPING 


In  order  to  complete  our  study  of  an  almost  symmetric  missile,  we 
will  now  assume  asymmetric  normal  force  and  damping  moments  in  addi- 
tion to  the  side  moment  we  omitted  in  section  III.  These  terms  are 
generalizations  of  those  of  Reference  9. 


CY  * 1 CZ  " 


- CN  6 - CN  C 
a a 


(7.1) 


C + i C = C tic  + [CMS  - i C M K 
m n m0  nQ  Ma 


+ i ^ c ♦ cM  u + cM  u 

a q q 


- i cM>  U'  * i a 

a 


+ 1 CM.  t«'  - 1 ♦' 

a 


(7.2) 


The  new  coefficients  are  defined  in  Table  IV.  This  assumed  force  and 
moment  can  now  be  inserted  in  Equations  (2.12-2.13)  and  y eliminated 
between  the  resulting  equations.  (Due  to  size  considerations,  products 
of  C*  terms  are  ignored.)  We  obtain 


MB  M f t'H'  - [Mr  + i (Mg  - <p'  H)]5 

= i A - H V - [M  + i (Mg  - H)]C  (7-3) 

where 

H ■ cn  - 2 CS  - K2  <c5  4 Si.) 

a q a 

MS  ■ P (CN  - CP  4 ki2  CMS 
a a 

" ■ * >C2  (Si  * s;  ) 


The  solution  is  assumed  to  have  the  form: 


(X,  + i u,)s  O,  + i “ )s 

S = k e 1 1 + k e 2 2 

l 2 


(X  - i w )s  Cx,  - i <*> Js 

+ k3  + k4  e l l + k5  e 2 ~ 


Substituting  in  Equation  (7.3),  we  have 


-i  { [Mr  - i (Mg  - ^ H)]  A - [M  + i (M  - <(/H)]a} 


M2  - M2  + (Mg  - H)2 


(Mg  - <p'  H)2 


kl  Cl  + 2 f 


p'  ^ - i (2  ♦'  - P)(Xj  - X1S)1 


k I M + i 


[(^  - *')H  + Mg] J = 0 


kj  j^M  + i [(Wl  + *')H  - Mg] J 

* \[C,  ♦ - P ■)(»,  - X1S) 

■ 2 Uj  (H  + 2 Xj  )]J  =0 


where 


- Mg] 

2<p'.  - P 


and  similar  equations  hold  for  k and  k . Outside  the  resonance 

2 5 

region,  the  real  part  of  the  determinate  of  the  coefficients  in 
Equations  (7. 6-7. 7)  set  equal  to  zero  gives  a slightly  modified  relation 
for  the  frequency  <f/  . The  imaginary  part  of  this  determinate  set  equal 

to  zero  gives  a good  approximation  for  the  damping.  With  the  approxima- 
tion 

C = - 2 f <t>'  w (7.8) 

l 1 


the  damping  exponent  has  the  following  simple  approximation  for  spins 
far  from  zero  and  the  resonance  region: 


X 

1 


X + 
IS 


M H 

f 4'  (2  ♦'  - P) 


(7.9) 


Similarly 


M H 

X = X _ + ; 

2 2S  f 4/  (2  4>  - P) 


(7.10) 


Thus  the  damping  exponent  differs  only  a little  from  the  damping  for  a 
symmetric  missile  with  the  average  moment  coefficients. 

The  bounds  of  the  resonance  region  are  set  by  the  zeros  of  the  de- 
nominator of  Equation  (7.5).  If  we  neglect  the  side  moment  coefficient, 
this  condition  has  the  simple  form. 


M2  - M2  + O')  2 F = 0 

where 

F = [H  - (2  - f) (C*  - C*)]2 

a 

- [H  - f C*  ]2 


(7.11) 


If  we  denote  the  damping-modified  resonance  spins  as 


bSR  and  ^aR 

and  assume  F to  be  small,  the  modified  resonance  spins  can  be  computed 

as  perturbations  of  the  zero-damping  resonance  spins  <)>7  and 

p 
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(7.12) 


r 


*.  ■[ 


4 (f  - i)  m re 


- r 

M 


4 (f  - 1) 


mJ  “ 


(7.13) 


Since  M was  selected  to  be  negative,  a positive  F reduces  the  size  of 
the  resonance  region  and  negative  F increases  it. 


VIII.  SUMMARY 

1.  The  general  motion  of  an  almost  symmetric  missile  is  well 
approximated  by  a symmetric  missile  with  average  coefficients. 

2.  Far  from  zero  spin  or  resonance  spin  rates,  the  first  observ- 
able modification  of  the  usual  tricyclic  motion  for  an  almost  symmetric 
missile  is  the  appearance  of  a 2<f>/  - (f/j  frequency  followed  by  the 

appearance  of  a 2<J>/  - <j/  frequency  as  the  asymmetry  becomes  greater. 

3.  Near  zero  spin,  both  of  these  additional  frequencies  have 
substantial  amplitudes  and  near  resonance,  the  24/  - 4/  frequency  has 

a substantial  amplitude. 

4.  For  spin  in  the  resonance  region,  large  trims  and  exponential 
undamping  are  possible. 


I 

f 
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LIST  OF  SYMBOLS 


k (C*  + i C*  ) 

t v m n 

0 0 

(f  <}/  )2/2  - M 


M2  - M2 
r 


Cf  Cm*  Cn 


drag  force 
P S V2/2 

missile-fixed  components  of 


^aerodynamic  moment 


p S £ Vz/2 


V CY’  Cz 


V C2 


C(  )’  C(  ) 


missile-fixed  components  of 


aerodynamic  force' 
, P S V2/2 


iii?  - f to.  + M , j = 1 , 2 
J J r i 


see  Table  IV  for  all  C definitions  not  given  above 


2 - <W 

-2S  [2  f ♦'  Mr  (“is  - “ns”’' 

»1S  P f <*.S  ' -is"'1 


-2  * 

CN  ” 2 CD  ' kt  °M  * CM.5 
a q o 


CN  + V <CM  + W 

a q a 


angular  momentum 

transverse  moment  of  inertia  ( = I = I ) 

y z 

axial  moment  of  inertia 


absolute  value  of  k. 

J 


absolute  value  of  k.. 

JR 
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LIST  OF  SYMBOLS  (Continued) 
complex  yaw  modes,  j = 1-5 

complex  yaw  modes  in  the  resonance  region,  j = 


V"*2 


reference  length 
k‘t2  CM 


k;2  CM 


M ♦ (f  - 1)(*')2 

p(c;  - * i;2  c* 

a a 


f *'  Cj 


<VV*' 


spin  rate  (rad/s) 

transverse  angular  velocity  components  in  the  missile- 
fixed  system 

reference  area 

dimensionless  arclength  along  the  trajectory 


u,  v,  w missile-fixed  components  of  V 

V missile  velocity 


X,  Y,  Z 


magnitude  of  V 

missile-fixed  axes  in  which  the  X-axis  is  along  the 
principal  axis  of  inertia  nearest  to  the  flight 
direction 
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X , X 
1 2 


LIST  OF  SYMBOLS  (Continued) 
tan'1  (w/u) , the  angle  of  attack 
sin'1  (v/V) , the  angle  of  sideslip 


/ 1 if  4j/  > <f/ 

(-1  if  0 < 4>/  < 4>'„ 


any  positive  number  small  in  comparison  with 

exponential  damping  coefficient  in  the  resonance-region 
solution  (5.1)  of  the  yaw  equation  (4.1) 

exponential  damping  coefficients  in  the  solution  (7.4) 
of  the  yaw  equation  (7.3) 

- (H<f/.  - M )/(2<f/.  - P),  j = 1,  2 
3 5 J 

(the  symmetric-missile  value  of  X^) 


(q  + ir)£/V 


v + 1 w 


, missile-fixed  complex  angle  of  attack 


£e1(*\  aeroballistic  complex  angle  of  attack 

air  density 

<p's,  the  roll  angle 

argument  of  k^R 

argument  of  k. 


plL/V  = spin  rate  (rad/cal) 

j-th  yaw  arm  frequency  in  the  aeroballistic  system, 
j = 1,  2 


(p  ± \/p2  4m)/2  , 

(the  symmetric-missi 


j = 1,  2 

le  value  of  <j>. ) 

J 


upper  boundary  of  the  resonance  spin  region 


LIST  OF  SYMBOLS  (Continued) 


u . 
3 


lower  boundary  of  the  resonance  spin  region 


j-th  yaw  arm  turning  rate  in  the  missile-fixed  system, 

j = 1,  2 

2 5* 

(_k  C*  ) , the  pitch  frequency  for  zero  spin 

ma 

-2  U 

(k^  C*  ) 2 , the  yaw  frequency  for  zero  spin 

8 

angular  velocity  of  the  missile-fixed  system  with 
respect  to  an  inertial  system 


Superscripts 

(') 

( )' 

( )* 
o 


= d ( )/dt 

= d ( )/ds  = (')A/V 
pSt 


2m 


( ) 


= complex  conjugate  of  ( ) 


t 

i 

! 
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APPENDIX  A.  EFFECT  OF  SMALL  MASS  ASYMMETRY 


In  this  report,  we  have  restricted  the  allowable  asymmetries  to 
the  aerodynamic  force  and  moment.  Mass  asymmetries  can  change  the 
center  of  mass  and  the  inertia  properties.  A transverse  shift  of  the 
center  of  mass  introduces  a trim  moment  due  to  drag  and  a roll  moment 
due  to  lift;  these  moments  can  be  easily  handled.1*  The  effects  of 
unequal  transverse  moments  of  inertia  and  of  products  of  inertia  are 
more  complicated  and  will  be  considered  here.  We  will  outline  the 
procedure  to  handle  small  inertia  asymmetries  and  show  their  impact 
on  the  results  of  this  report. 

For  asymmetric  inertia  properties,  the  angular  momentum  vector 
has  missile-fixed  components 


H = I p - J q-J  r 
x x r xy  1 xz 


(Al) 


H = - J p + I . [ Cl  + d.)  q + dr] 
y xyr  t 1 ln  2 


(A2) 


V'.:'*  't  td2«  * (1  - V r> 


(A3) 


where 


dl  ’ <Jy  - V/(2  V 


d.  = - J /I 
2 yz  t 


*t  = 


12.  D.A.  Price,  Jr.,  " Sources , Mechanisms,  and  Control  of  Roll  Resonance 
Phenomena  for  Sounding  Rockets, " Journal  of  Spacecraft  and 
Rockets  £,  November  1967,  pp.  1516-1625. 
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and 


J , J , J are  products  of  inertia, 
xy  / z xz 


Equation  (2.13)  now  takes  on  the  form 


U7  + i (4/  - P)u  + (d1  + i d2)  (tj7  + i <p'  y) 

= k“2  (C*  + i C*) 
t m n' 


(A4) 


+ C*  [u  + (di  + i d2)p] 


* 1 ‘J,y  * 1 


Cp  can  now  be  neglected  and  the  static  moment  of  Equation  (3.5) 
can  be  inserted  to  obtain  a new  form  of  Equation  (4.1): 


+ i f <t>'  £'  - M K + M l 
r 


(A5) 


- (dx  ♦ i d2)[l"  * («')2C]  = i 


where 


A = k'2  (C*  + i C*  ) + i (j  + i J )(<j/)2/l 

t m n xy  xz-'  J 


If  the  pentacycle  solution.  Equation  (4.4),  is  substituted  in 
Equation  (A5) , we  obtain  the  following  relations  for  the  modal  ampli- 
tudes : 


k3  = - 1 [Mr  A - [M  ~ (di  + 1 d2)  (♦'  )2]Ajb_1 


(A6) 
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= [m  ♦ (d1  ♦ i d2)[u2  - (^)2JJkj  c;1 


= + (dj  + i 


i d2)[u2  - (*')2]  k2  C2'J 


where 


C?  + 2 f u.  C - £m  + dx  [w?  - (<t/)2]J' 


* d2  [oj?  - (4/  )2]  =0 


C.  = a»?  — f <J>/  to . + M 
J J J r 


b = Mj  - [M  - di  (^)2]Z  - d2  (4/)  4 


It  should  be  noted  that  for  general  motion  away  from  zero  spin  or  the 
resonance  region,  a good  approximation  for  is  still  - 2 f ^ u.  . 

The  new  frequency  equation  can  be  obtained  from  Equation  (A9)  by 
eliminating  C.  and  assuming  small  mass  asymmetry  (i.e.  neglecting 

d2  in  comparison  with  one) . 

u)4-2aio2  + b=  0 (A10) 

J J 


where 


a « (f$')2/2  - Mr  + M dj  - (d2  + d2)(4>'): 


The  pitch  and  yaw  frequencies  for  zero  spin  are  given  by 
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u2  = - k'2  C* 
a y m 

J a 

(All) 

w2  = k-2  C* 

6 z n8 

(A12) 

where 

k = /i  /m  £2 

y y 

k = /i  /m  a2 
z z 

The  coefficients  of  the  frequency  equation  can  be  expressed  in  terms 
of  these  zero-spin  frequencies. 
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